Jahn- Teller systems from a cavity QED perspective 
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Jahn- Teller systems and the Jahn- Teller effect are discussed in terms of cavity QED models. By 
expressing the field modes in a quadrature representation, it is shown that certain setups of a two- 
level system interacting with a bimodal cavity are described by the Jahn- Teller Exe Hamiltonian. 
We identify the corresponding adiabatic potential surfaces and the conical intersection. The effects 
of a non-zero geometrical Berry phase, governed by encircling the conical intersection, are studied 
in detail both theoretically and numerically. The numerical analysis is carried out by applying a 
wave packet propagation method, more commonly used in molecul or chemical physics, and analytic 
expressions for the characteristic time scales are presented. It is found that the collapse-revival 
structure is greatly influenced by the geometrical phase and as a consequence, the field intensities 
contain direct information about this phase. We also mention the link between the Jahn- Teller effect 
and the Dicke phase transition in cavity QED. 

PACS numbers: 42.50.Pq, 03.65.Vf, 31.50.Gh, 71.70.Ej 



I. INTRODUCTION 

The Jahn-Teller (JT) effect, due to Hermann Jahn and 
Edward Teller states that a symmetry-breaking is 
likely (only exceptions are linear molecules or molecules 
possessing Kramers degeneracy points [2] ) to occur when- 
ever there is an isolated degeneracy of electronic states 
in a molecule, a so called conical intersection (CI) 
Over the years, the JT effect has gained enormous atten- 
tion, mainly in molecular and condensed matter physics 
@, IE @] ■ A simple model system Hamiltonian possess- 
ing a CI, later termed Exe, was presented by H. C. 
Longuet-Higgens et al. The main result of this work 
states that the angular momentum quantum number is 
half integer valued rather than an integer. This phe- 
nomenon arises from a geometric phase, on top of the 
dynamical one, obtained while encircling the CI. The 
additional phase must be introduced in order to have 
a single valued total (electronic and vibrational) wave 
function. This was further analyzed in Q, where it was 
shown that the double valuedness of the electronic wave 
function can indeed be removed by introducing a "vec- 
tor potential" term in the Hamiltonian. For CI models, 
this resembles the Aharanov-Bohm effect 9], and gave 
rise to the molecular Aharanov-Bohm effect and molecu- 
lar gauge theory [ljl HH • A deeper understanding of this 
phase effect was gained with the seminal paper of M. V. 
Berry [12j | , which presents a general formalism for the ge- 
ometrical phase factors that an adiabatic change in the 
Hamiltonian brings about. In the spring from this work 
came several papers on the geometrical phase related to 
CIs, see Refs. (flEI- 

The effect of the geometrical phase on physical ob- 
servables has since then been discussed and experimen- 
tally verified in several reports [HI, [3| ■ The modulation, 
caused by the geometrical phase, of the wave function 
has been addressed in Ref. [1 51 ] , where a dynamical wave 
packet approach, like the one used in this article, is ap- 
plied to the Exe JT model. Inclucion of spin-orbit cou- 



plings has been investigated in terms of the JT effect [16 



and of wave packets [lTf . Recently, other properties of 



JT models, not only the Exe, have been considered, for 
example, quantum chaos (l9| and ground state entangle- 
ment 

Although the JT effect has not, to the best of my 
knowledge, been discussed within cavity quantum electro- 
dynamics (QED), the geometrical Berry phase has been 
analyzed in the framework of cavity QED [U [52, [HI- 
These references study the effects induced by the vac- 
uum field on the geometrical phase. In other words, the 
treatment of the two-level particle in the time-varying 
field is then considered on a fully quantum mechanical 
footing. The degeneracy point is not, however, identified 
as a CI in these works, and the rotating wave approxi- 
mation (RWA) has been applied which is not the case in 
the present article. In addition, in the situation studied 
here, the geometrical phase is said to be of dynamical 
character as it originates from the intrinsic evolution of 
the system rather than from "external" changing of 
the Hamiltonian [U [2^, [23j] . Thus, the circumstances 
and approaches are not ably different between this work 
and the ones of Refs. [U [22], [23|. Among others, we 
examine the cavity QED system in a conjugate repre- 
sentation in which the intracavity fields are expressed in 
their quadrature operators rather than the standard used 
creation and annihilation ladder operators. In this pic- 
ture, the link to JT systems and to CIs is revealed and, 
in fact, the Dicke normal- superradiant phase transition 
in cavity QED is seen to be related to the JT effect. 

The model system is a two-level quantum-dot embed- 
ded in a cavity and interacting with two degenerate field 
modes. The numerical analysis is carried out using wave 
packet propagation methods; an initial state of the sys- 
tem is let to evolve under the corresponding Hamilto- 
nian. The full Hamiltonian dynamics is considered and 
compared to results obtained from a second Hamiltonian 
which shares the same adibatic potential surfaces (APS) 
but lacks a geometrical phase. We give analytical ex- 
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pressions for the characteristic time scales, collapse and 
revivals, for both systems and it is found that the revival 
time is in a sense twice as long in the case where the 
geometrical phase is excluded. A consequence of this is 
reflected in the intracavity field intensities. Thus, mea- 
surements of the field intensities of the two modes give a 
direct indication of the geometrical phase. 

The outline of the paper is as follows. In Sec. |TT] we 
review some basics of JT models and CIs. Especially, in 
Subsec. Ill Al we introduce the Exe Hamiltonian and dis- 
cuss its APSs, while Subsec. Hi Bl derives the geometrical 
phase accumulated by encircling the CI, and Subsec. Hi CI 
considers the JT effect in general and in the Exe Hamil- 
tonian in particular. The next Sec. IIIII is devoted to our 
cavity QED model and it is shown how the Exe model 
occurs for a two-level system interacting with degenerate 
bimodal fields. A discussion of the corresponding JT ef- 
fect in cavity QED is outlined in Subsec. IIIIB| where a 
parallel with the Dicke normal-superradiant phase tran- 
sition is drawn. Our numerical results of the cavity JT 
system are presented in Sec. lIVI both analyzing the short 
and long term behavior and how the geometrical phase 
comes into play. Finally we summarize in Sec. fVl 



II. THE JAHN-TELLER MODEL 

J ahn- Teller systems are characterized by a degeneracy 
point of coupled potential surfaces, a CI. In one dimen- 
sion, the simplest example is the Ex/3 model, also called 
Rabi or spin-boson model [24j |. It describs a spin 1/2 
particle coupled to a single boson mode [25[. In certain 
parameter regimes, a RWA can be applied in which this 
model relaxes to the one of Jaynes and Cummings [2t| . 
In one dimension, the wave packet (state of the system) 
cannot encircle the CI without passing through it, and 
therefore there is no corresponding dynamical geometric 
phase [12|. In two dimensions, generalization of the Ex (3 
model leads in certain situations to the Exe model which 
will be the subject of this section. 



A. The Exe Hamiltonian 

The simplest Jahn- Teller Hamiltonian with two vibra- 
tional degrees of freedom is the so called Exe one, given 

by 



H 



JT 



2m 2 



2m 



(1) 



+Xxa x + Xy&y. 



Here pi and Xi are momentum and position in the i di- 
rection of the "particle" with mass m. The <r-operators 
are the standard Pauli matrices obeying the commuta- 
tion relations 



|<T,-,<7j 



(2) 



and with the z-eigenstates, 6" z |±) = ±|±) and A the cou- 
pling constant. Clearly, at the origin x = y = 0, the two 
potential surfaces are degenerate. In the presence of ei- 
ther spin-orbit coupling [TtI HtJ or an external magnetic 
field (20l.[28j. an additional detuning term Aa z /2 is added 
to the Hamiltonian, where A is the spin-orbit splitting or 
the magnetic strength. With this term present, the de- 
generacy is lifted and the intersection becomes avoided. 

The form of the Hamiltonian (jTJ) defines the diabatic 
basis and diabatic potentials, namely; a diabatic state is 
written as ^(x,y) — f±(x,y)\±) for some normalized 
function f±(x,y), and the diabatic potentials, once the 
detuning A is included, are muj 2 (x 2 +y 2 ) /2±A/2. Before 
defining the APSs, we express the Hamiltonian in polar 
coordinates 



x ± iy = pe ±itp 



(3) 
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Id Id 2 



2m \dp 2 p dp p 2 dtp- 2 



A 
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Xpe^ 
A 
~ ~2 



(4) 



Let us introduce the unitary operator [8|, [if 



U 



sin(V) cos(V) 
cos(^)e IM sin(V)e l/1 



where 



tan(2f) 



2Xp 



(5) 



(6) 



and ji = ip, which diagonalizes the last term of Eq. ([?]). 
However, U does not commute with the kinetic term in 
Eq. (j4|) and consequently, the transformed Hamiltonian, 
Hjt = U~ 1 HjtU , is non-diagonal. The off-diagonal 
terms are the non-adiabatic couplings, which usually are 
small far from the crossing. Omitting these terms defines 
the adiabatic Hamiltonian 



H% = T + VI 



ad 



v ct 



V 



gauge 3 



(7) 



where 



T 



11 



d 2 



Id Id 2 



2m \dp 2 pdp p 2 dip 2 




(8) 



Yr, 



2{& 2 +AX 2 p 2 Y 



V gauge — ^ 2 ' ^ 



A 



2p 2 \ ^A 2 +4X 2 p 2 A 2 dip 



1 d 
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We will assume that the evolution takes place mainly 
upon the lower APS, i.e., a z = —1, but in the numer- 
ical simulations we consider full dynamics without any 
approximations. Close to the crossing, the non-adiabatic 
couplings may have a significant impact on the dynamics 
[Jill]. The term V£ d defines the APSs, while V cen t and 
V gauge, arising from the commutator between U and the 
kinetic energy operator, are centrifugal corrections [30| . 
The last term has been labeled gauge for reasons that will 
become clear later on. We display two examples, A = 
and A ^ 0, of the APSs V± d in Fig. [TJ For a non-zero 
detuning, as pointed out, the crossing at the CI becomes 
avoided, with splitting amplitude A. The lower surface 
has the familiar sombrero shape, while the upper pos- 
sesses a single global minimum at x — y = 0. For large 
detunings A, the Mexican hat structure is lost, and the 
minimum of the lower APS is at the origin. Especially, 
the radius giving the potential minima is given by [17| 




FIG. 1: Two examples of the adiabatic surfaces in the Exe 
Jahn- Teller model with detuning A = (a) and (b). 



B. Berrys geometrical phase in the Exe model 

The adiabatic states, defined by U, are arbitrary up to 
an overall phase. The phase choice in Eq. ((5]) is cho- 
sen such that the states are singled valued as ip is varied 
by 2tt. For example, the alternative obtained by multi- 
plying U by exp (—iip/2) implies double- valued adiabatic 
states. Unitary transformation of the Hamiltonian in this 
second case leads to an adiabatic Hamiltonian lacking the 
term proportional to id /dip. The very last term of iff]), 
containing the differential operator d/dtp, can be viewed 
as a vector potential. Indeed, this term can be combined 
with the canonical momenta to define a kinetic momenta. 
Thus, the two options of overall phase of the adiabatic 
states given above result in either single-valued states 
with a vector potential present in the Hamiltonian or no 
such vector potential term but double- valued states [3l[ . 
The source of a vector potential term is the cause for 



having a non-zero geometric phase as the system encir- 
cles the CI in analogous to the Aharanov-Bohm effect 

For the system evolving along a closed loop C in pa- 
rameter space, the geometrical phase can be calculated 
according to QJJ 

7 „(C) = £<n(R)|V R n(R)) • dR. (10) 

Here, |n(R)) is the n'th adiabatic eigenstate and R 
the set of parameters. In particular, in our case 
|n(R)) = (sm(9), — cos(9)e lv ) and as we consider a time- 
independent dynamical problem, the varying parameters 
R are the coordinates p and ip. Especially, we consider a 
wave packet located at the minima of the sombrero po- 
tential, such that (p) w Pmin, and ip is changed from 
to 2ir. For a general radius R we find the Exe geometric 
phase (l7| 

Ut{R) = -7T [ 1 + 1, (11) 

which for R = p m i n becomes 

iMPmin) = -7T (l + ^ V (12) 

For A = we obtain the well known sign change of the 
wave function when encircling a CI, causing half integer 
angular momentum quantum numbers. 

C. The Jahn- Teller effect 

Using group theoretical arguments, Jahn and Teller 
proved that for almost any degeneracy (CIs) among elec- 
tronic states in a molecule, a symmetry-breaking is "al- 
lowed" which removes the degeneracy and lowers the to- 
tal energy of the model system ground state [l[ . It turns 
out that this symmetry-breaking indeed takes place in 
the majority of cases with some exceptions [U, Hence, 
the molecule favors a distortion of its most symmetrical 
state. This effect is quenched when spin-orbit coupling 
is taken into account, but may still exist [lq |. Returning 
to the APSs of the Exe Hamiltonian ([7]), the state with 
highest symmetrical is a wave packet centered around the 
origin (x = y = 0). Loosely speaking, for vanishing de- 
tuning A, is it intuitive to expect the wave packet to slide 
down the potential surfaces towards the minima of the 
sombrero. Semi-classically, the wave packet experiences 
a non-zero force F = —W(x,y). For large detunings 
however, we saw from Eq. @ that the potential surfaces 
may possess a single global minimum which will prevent 
the JT-distortion. On the other hand, for small but non- 
zero detuning, the sombrero structure is present for the 
lower APS and here quantum fluctuations will permit a 
symmetry breaking. Naturally, the above arguments are 
semi-classical and the full evolution is quantum mechan- 
ical and described by the coupled system. Nonetheless, 
it gives some insight and intuition of the dynamics. 
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III. JAHN-TELLER MODELS IN CAVITY QED 

Spin boson models naturally occur in cavity QED. 
Here the boson subsystem represents a single or several 
quantized modes of an intracavity field, while the spin de- 
grees of freedom describes either two- level atom/atoms 
[32| or solid state quantum-dot/dots 33]. Contrary to 
standard formulations of Jahn- Teller models, here the 
bosons are the photons of the field rather than vibrational 
phonons, and the internal structure corresponds to two 
discrete energy levels of the atom or quantum-dot. In the 
single mode case, a microscopic derivation gives the Ex (3 
Hamiltonian [35j ] in the assumption of dipole approxima- 
tion and neglecting the self-energy (see below). In most 
cavity QED experiments involving atoms, the applica- 
tion of the RWA is justified, in which the Hamiltonian 
identifies the analytically solvable Jaynes-Cummings one 
[26| . The APSs (or rather adiabatic potential curves) of 
the Jaynes-Cummings model contain the differential mo- 
mentum operator [2J|, and they are therefore said to be 
of non-potential form. Nonetheless, even if the picture 
of potential surfaces is less intuitive due to the momen- 
tum dependence, the JC model renders a sort of general- 
ized CI (curve crossing). To go beyond the RWA regime, 
the coupling to the field must be substantially increased 
compared to atomic cavity QED setups. This is indeed 
the case for solid state quantum-dots coupled to a cav- 
ity. In fact, the crucial parameter, coupling divided by 
the two-level transition frequency, can be made several 
orders of magnitude larger in condensed matter systems 
[H| compared to atom-cavity ones. Another possibility 
to achieve ultrastrong atom-field couplings is to consider 
Bose-Einstein condensates coupled to an intracavity field 

tu. 



A. The model system Hamiltonian 

To obtain CIs rather than curve crossings, multi-mode 
cavities must be considered 36]. For simplicity we will 
assume two degenerate cavity modes such that they share 
the same frequency uj and also same coupling amplitude 
A to the quantum-dot. The Hamiltonian in the dipole 
approximation reads [3?], HH 



H cav = IvjJ {a^a + Pbj + fo-^&z 



(fit + a) (a+e-** + as**) ( 13 ) 



+ (fit + b) {a+e- w + a~e w ) 



Here at and b' (a and b) are creation (annihilation) oper- 
ators for the two field modes, (f> and 9 field phases, f2 the 
quantum-dot transition frequency, and 2<r ± = a x ± ia y . 
In the following we will label the two cavity modes by 



a and b. Before proceeding, for brevity we introduce a 
characteristic energy hu> and time scale such that 

we consider dimensionless variables 



A 



A 



n 



r = tit, 



(14) 



where t is the unsealed time. In a conjugate variable 
representation defined by the operator relations 



1 



Px 



Py 



ftt 



1 



V2 



7^ + h \ 



(15) 



v/2 



where [i^pj 
the form 



[y,p y ] = i, the Hamiltonian (|13p takes 



fy 



X 

T 



2 



2 2 2 2 2 
+2Ax[cos(<^>)<7 x + sin(0)<T y ] 

+2A$ [ cos(#)a x + sin(0)<7 a ] . 



(16) 



For the simple example of <j> = and 9 — ir/2 we recover 
the cylindically symmetric Exe Hamiltonian (JTJ) . In fact, 
for 



|^-0|=O' + 1/2)7T, 



j integer, (17) 



Hcav is unitarilly equivalent with the Exe Hamiltonian 
Hjt by identifying with A. In some special cases of 
the phases, the last two terms of (jTHJ) can be written as 



2X(x + y) 



2X(x - y) 



e _< * 
e^ 











for 9 



for 9 



2j7T, 



(2j + 1)7T, 



. . (18) 

for some integer j. For these situations the CI is re- 
placed by an intersecting curve in the directions of ip = 
37r/4, 77r/4 or cp — tt/A, 5n/A respectively. Here, it is 
clear that by a unitary rotation, the adiabatic states can 
be made real, indicating that the geometrical phase be- 
comes identically zero as the wave packet is encircling the 
CI. This is indeed seen in Fig. [5] displaying the geometric 
phases of the Hamiltonian (|16p . The general form of the 
APSs, in polar coordinates, reads 




4AV[i + cos(0-(9) sin(2(^)]. 



(19) 

The lower APS has two minima for angels ip = 7r/4, 37r/4. 
It is known that the "quadratic" Exe Hamiltonian has 
three local minima in the sombrero shaped potential [l5l . 
[13] ■ This derives from a term of the form sin(3c/?) in the 
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APSs. Here we have a sin(2(/?)-dependence instead and 
hence the double minima structure. 

The single valued adiabatic states can again be written 



*«(P, <P) 



but with 



sin(j/) 
cos(^)e l/i 



4Xp 



cos(i/) 
sin(^)e l/i 



tan(2i/) = y/l + sin(2^) cos(0 - 9), 



H — tan 



cos(y>) sin(0) + sin(y>) sin(0) 
cos(<p) cos(0) + sin(y>) cos(#) 



(20) 



(21) 



Encircling the CI at a radius i? we find the geometric 
phase (fH 



7cow(-R) = - 



cos (a)-^— ctw 



(22) 



p=R 



The phase (|2"2"j) is depicted in Fig.[5]as a function of 9 and 
A for fixed — and Q = 1 (a) and = 1/2 and S7 = 1 
(b). The asymptotic value for large couplings A is either 
— 7r or (modulo 2ir). The radius R is taken to be the 
minimum of the cylindrically symmetric case (9 = 7r/2). 
For = 0, the greatest effect generated by the geometri- 
cal phase is seen to be in the symmetric case of 9 — n/2, 
while for 0^0 the situation becomes more complex. 
Note that, according to Eq. ([9]), and identifying f2 with 
A, A ^ in order to have a sombrero structure of the 
lower APS, which is the reason why A is not approaching 
0. 



B. The Jahn- Teller effect in cavity QED 

The Jahn- Teller effect states that, in the presence of 
a CI, lowering the symmetry may be energetically fa- 
vorable in various systems, typically for molecules and 
crystals. In terms of the Exe Hamiltonian this is intu- 
itive, since a wave packet centered at the origin (the CI) 
will have a larger "potential energy" than a wave packet 
located at p m in- The corresponding symmetry breaking 
in cavity QED implies that the system ground state con- 
sists of non-zero fields in the two modes. At the same 
time, the quantum-dot is not entirely in its lower state 
but in a superposition of its two internal states. This is 
related to a well known phenomenon in quantum optics, 
namely superradiance [391 ]. This, of course, comes about 
due to the strong interaction between the quantum-dot 
and the cavity fields. For multi quantum-dot systems 
and in the thermodynamic limit, where the number of 
two-level quantum-dots and the volume tend to infinity 
while the density is kept fixed, this results in a second 
order quantum phase transition between a normal (the 
field in its vacuum) and a superradiant phase (a macro- 
scopic non-zero field) |4(j. The critical coupling of this 
phase transition is given by A c = yj\ A\uj 40], which in- 
deed follows from Eq. ©. 





FIG. 2: The geometric phase (|22[) as function of coupling 
A and field phase 6. The radius R = pmin, where p m in = 
\/4A 2 — (fi/4A) 2 is the adiabatic potential minima of the som- 
brero when 9 — <j> — tt/2 (the cylindrically symmetric case). 
The other dimensionless parameters are = and £1 = 1 (a) 
and = 1/2 and fi = 1 (b). 



However, it can be shown that for a quantum-dot in 
which the lower state is its ground state, the normal- 
superradiant phase transition is an artifact from neglect- 
ing the self-energy term from the Hamiltonian [4l| . For 
a single mode, and in unsealed units, this term is given 
by 



e nh 
2m Via 



(23) 



where e is the electron charge, m its mass and V the 
effective mode volume. This should be compared with 
the matter-field coupling 



x = nd 



2-Kh 



(24) 



where d is the dipole moment of the transition of inter- 
est in the quantum-dot. It is clear that the self-energy 
term H se tends to quench the Mexican hat structure, 
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and further that A and H sc are not fully independent. 
Indeed, in Ref. [4l[ it is demonstrated, either using the 
Thomas- Reiche-Kuhn sum rule or simple thermodynam- 
ical and gauge invariance arguments, that the sombrero 
shape cannot be obtained for any set of physical param- 
eters. 



To circumvent this obstacle one may use a two-photon 
Raman type of interaction, where three levels of the 
quantum-dot are coupled through the cavity mode and 
an external classical laser field 0, EH . In the large de- 
tuning limit of one of the internal levels it can be adia- 
batically eliminated [44| and one arrives at an effective 
model very similar to the one above. In such procedure 
one introduces an additional independent parameter, the 
detuning S of the eliminated level [45[ . As the detuning 
enters in the effective matter-field coupling parameter, 
but not in the self energy term, these two become in prin- 
ciple independent, and in particular A can be made large 
in comparison with H se . The effective model, once the 
detuned level has been eliminated, contains some Stark 
shift terms that will modify the potential surfaces, but 
the CI and sombrero structure are still present. The 
external laser fields have the advantage of being easily 
controllable in terms of system parameters such as am- 
plitude and phase. A system Hamiltonian suitable for re- 
alizing the Jahn- Teller model can be found in Ref. [22| . 
It should be noted though, that the effective Hamilto- 
nian is in general time dependent, which is prevailed by 
imposing a RWA. This, however, induces a "momentum- 
dependent" potential surface, but nonetheless, the JT 
symmetry breaking is still present in this approximation 

m 



Another possibility to surmount the problem with the 
self-energy, which is assumed in this paper, is to use 
the fact that the two internal levels that couple to the 
cavity modes are normally highly excited meta stable 
(Rydberg) states. For these states, neither the Thomas- 
Reiche-Kuhn sum rule nor the thermodynamical argu- 
ments apply, and the symmetry breaking may still oc- 
cur. We therefore discard the self-energy terms as they 
would only modify the frequencies of the harmonic po- 
tentials. In general, also for the Raman coupled model, 
the states involved are highly excited metastable states 
and the sum rule cannot be applied in those cases ei- 
ther. A benefit of the Raman model, compared to a one 
photon model, is the higher controllability of the system 
parameters; especially the diagonal element is a detuning 
parameter (and not a transition frequency f2) that can be 
made small compared to the matter-field coupling. The 
drawback of an effective Raman model is that the anal- 
ysis is considerably less intuitive due to the RWA. Here, 
we therefore choose the simpler model of the two as the 
physical phenomena may be more easily extracted from 
it. 



IV. NUMERICAL RESULTS 

Contrary to molecular or solid state systems, proper- 
ties of the cavity fields are in comparison easily measured, 
for example, phoson distribution [4^|, field quadratures 
[47l | and, in fact, the whole phase space distribution us- 
ing quantum tomography [48|. Using wave packet prop- 
agation methods, the dynamics of such quantities will 
be studied in this section with emphasizes on the effects 
emerging from the geometrical phase. As an initial state 
we take a disentangled one, given in cartesian coordinates 
by 



9(x,y,0) =tp(x,y,0) 



1 

71 



1 

-1 



(25) 



where 



\/7T 



-(lmxa) — (Imj/o) 



(26) 

The initial quantum-dot is a linear combination of its two 
internal states with equal amplitudes, and the two field 
modes are in Gaussian states corresponding to coherent 
field states; |afo/v2) and \yo/V%) respectively [50]. Such 
initial states are readily prepared experimentally. We 
will further pick yo = and xq = 2 A such that the initial 
wave packet is approximately centered at the minima of 
the sombrero. Note that the initial average momentum 
is zero, and that ^(x, y, 0) is different from the adiabatic 
states ([20)) . A consequence of this is that the wave packet 
evolution will not be restricted to a single APS. However, 
the upper adiabatic state is only marginally populated 
for our particular choice of initial state and the main 
phenomena studied here, the effects of the geometrical 
phase on the field properties, is indeed seen even though 
slight interference between the two adiabatic states oc- 
curs. Hence, we emphasize that the dynamics take place 
mainly on the lower adiabatic surface. Properties of the 
upper APS have been studied in Ref. [49(. 

We restrict the analysis to the cylindrically symmetric 
case, where the time evolved state 



$(x,y,T) 



1 

71 



will predominantly spread along the minima of the som- 
brero potential. As the wave packet broadens it will, after 
a certain time, start to self-interfere. We may estimate 
the characteristic time for this process by approximate 
the inherent spreading by free evolution along the min- 
ima of the sombrero potential to get 



T 



V4tt 2 A 2 



1 ps 2ttA. 



(28) 



Within this time, the wave packet width has expanded 
over a distance 2irp m i n . 

From the full system state (|2"T|) . we can derive the re- 
duced density operators for the separate constitutes 



Pi(r) = Tr i; fe p(r) 



(29) 
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where the subscripts represent, either the 
two modes a and b or the quantum-dot, and 
p(r) — ^*(x,y,T)^(x,y,T), Using the reduced 
density operators we will especially study the photon 
statistics and the Husimi Q-distribution [111 HH 

Pi(n) = (n\pi(r)\n), 
(hi) =} nPAn), 

V ( 3 °) 

Qi{a) = —(a\pi(T)\a). 

7T 

Here, \n) is the n'th-photon Fock state, \a) a coherent 
state with amplitude a and the subscript i = a,b for the 
respective modes. 

A. Dynamics on the Ti n time scale 

Discussed in Sec. [TlJ it is the term V gauge that gives 
rise to a geometric phase. To correctly describe the adi- 
abatic evolution each term of -Hjf- in Eq. ([7]) must be 
taken into account, and it is not enough to study dy- 
namics upon the potentials V± d . In this subsection we 
study the full dynamics using Hamiltonian (fTB]) . and we 
hence go beyond any adiabatic approximation. However, 
in order to identify the effects of the geometrical phase 
we compare the results with the ones obtained by prop- 
agating the same initial state using the "semi-adiabatic" 
Hamiltonian defined as 

H# = T + (31) 

where T and V^ d are both given in Eq. ([5]) . Accordingly, 
a wave packet evolving via the Hamiltonian Hj? around 
the origin will not accumulate any geometrical phase. 

The characteristic time scale Ti n determines how long 
it takes for the particular initial state PS]) to inherently 
spread out across the CI and start to self-interfere. It is 
therefore a measure of the collapse time. The effect of 
the geometrical phase on the probability wave functions 
|*(x,?/,t)| 2 = \Mx-JLT)^/2 + \^ g (x,y,T)\ 2 /2 has been 
discussed in Refs. [15|, QjJ. Initially we choose xo > 
while j/o — such that interference of the evolved wave 
packet sets off at —xo where the two tails of the packet 
first join. Destructive and constructive interference cause 
nodes (vanishing probability distribution) and anti-nodes 
(non-vanishing probability distribution) in |\I/(x, y, t)\ 2 , 
and the ring-shaped wave packet splits up in localized 
blobs. In the case of Hjt, in which the geometric phase 
is zero, an anti-node builds up at x = — xq, while for 
^jriPmin) = (as is the case of zero detuning in the 
Exe model) a node is formed at x = — xq. The location 
of the corresponding node or anti-node depends on the 
value of jjT(Pmin), and in all our examples f2 < A such 
that "fjT{pmin) ~ giving a node at x sa — xo- These 
features are visible in Fig.[3]showmg the numerical results 



of the propagated distributions \^(x, y, r)| 2 for three dif- 
ferent times r. Full dynamics governed by Hamiltonian 
(fT5|) . with </> = and 6 = ir/2 (cylindrical^/ symmetric 
case) are shown in the left plots, while the right ones 
reproduce the results from propagation using the semi- 
adiabatic Hamiltonian (|3 1 1) . The effect of the geometri- 
cal phase becomes clear once the wave packet starts to 
self-interfere. The number of localized blobs depends on 
time r and system parameters and then especially p m in- 
Note that very similar results where presented in Refs. 
[HI However, the Hamiltonians used for the sim- 

ulations in Refs. [lol QjJ are in general different from 
the one H°A, utalized here; a single surface approxima- 
tion [3, [52J is applied in most examples of Refs. [H, [TtJ 
while here it is only considered for the non-geometrical 
phase case. 



With geometric phase Without geometric phase 




-5 5 -5 5 

x x 

FIG. 3: Snapshots of the wave packet distributions 
|*(x,2/,t)| 2 at times r = T in /4 (a) and (b), r = T m /2 (c) 
and (d), and r = Ti n (e) and (f) for the cases with (left) and 
without geometrical phase (right). In the last two plots, the 
difference between the interference structures is clearly visi- 
ble. See Refs. [l5l . Il7j for similar results. The dimensionless 
parameters are SI = 0.5 and A = 3. 



As the initial wave packet starts to spread, a non-zero 
field will begin to build up in the vacuum b mode, on 
the cost of decreasing field intensity of mode a. How- 
ever, without the RWA, the total number of excitations 
is not conserved. In Fig. 3] we display the individual pho- 
ton distributions Pi(n) at a quarter of the interference 
time Ti n . Already at this instant has the initially empty 
mode a non-zero field intensity, and its photon distri- 
bution consists mostly of even number of photon states. 
This is a typical characteristic of Schrodinger cat states 
[32| , and in the next subsection we will indeed show that 
such a state is created in the system at certain times. 
The small but non-zero population of odd photon num- 
bers in the b mode is caused by non-adiabaticity; for the 
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semi-adiabatic Hamiltonian (|3ip the odd photon num- 
bers are never populated for the given initial state ([25]) 
with yo = 0. 



0.5 




5 10 15 20 25 30 

n 

FIG. 4: The photon distributions P a (n) (white) and Pb(n) 
(black) after a time r = Tin/4. The dimensionless parameters 
are the same as in Fig. [3] 



B. Dynamics beyond the T„ time scale 

Seen in the previous subsection, for an initial local- 
ized state mainly located at the minima of the som- 
brero potential and with zero average momentum, the 
time scale Tj„ determines the collapse time; the time it 
takes for the localized wave packet to spread out over 
its accessible phase space region. Over longer periods, a 
revival structure in physical quantities is expected [53j, 
where localized bumps are formed in phase space signal- 
izing fractional or full revivals (54j. It has been pointed 
out however, that the collapse-revival characteristics are 
rather different in models where the RWA has been ap- 
plied [H, [EH . Typically, in the RWA regime the various 
time scales become long. In this work we are outside such 
a regime, and we will in particular find that the revival 
time is given by a multiple of ATi„ and that phase space 
evolution is significantly different for the two Hamiltoni- 
ans ([To]) and l|3"Tj) due to the geometrical phase. 

From Fig. [3] we see that after a time Tj n , the initial 
wave packet is spread throughout the minima of the som- 
brero potential and the self-interference causes nodes in 
the probability distribution. The number of localized 
bumps depends on p m in ©, but also on the time r; 
at first, when the self-interference sets off, the number 
of bumps increases to a maximum value and then the 
number begins to decrease and eventually form a single 
localized wave packet. Full revival occurs when a single 
localized bump is formed at the same position as the ini- 
tial wave packet. To study the field dynamics we use the 
Q-function for the two modes a and b of Eq. ([50)) . We 
will present the two functions Q a and Qb in the same 
plots for brevity, but mark them with letters a and b 
respectively. At the initial time t = 0, Q a and Qb are 
Gaussians centered at a = (y/2X, 0) and a = (0, 0) re- 
spectively. As time evolves, the b mode builds up its 



intensity and the Q-function moves away from the ori- 
gin, while Q a at first decreases its intensity by tending 
towards the origin. However, over longer time scales, 
t > Ti n , a swapping of energy between the two modes 
will take place. This phenomenon has been discussed 
in our model, but only when the RWA has been imposed 
[HI . As our analysis concerns a regime far from the RWA 
one, this exchange of energy between the modes occurs 
at very different time scales than in Ref. [38j], similar to 
what was found for the inversion in the JC model [24| . 
Namely, the characteristic time scales in the parameter 
regimes of the RWA and without the RWA in the JC 
model can differ by orders of magnitude. 

From our numerical simulations we have found that lo- 
calization of the phase space distributions comes about at 
multiples of time Tf rac — ATi„ , which hence are the char- 
acteristic scales for fractional revivals [EH, . The larger 
the radius p m in, the better resolved wave packet localiza- 
tions. In Fig. [5] we display examples of the Q-functions 
Q a and Qb (indicated in the figures by a and 6) obtained 
either from the full system Hamiltonian (|I6[) (left) or 
from the semi-adiabatic Hamiltonian (f3Tj) (right). The 
times are here, r = ATj„, 2ATi„, 3ATi„, 4ATi„. A clear 
discrepancy is seen between the two models. For exam- 
ple, at t = 2XTi n (c) and (d), mode a in the left plot (with 
geometrical phase) is approximately in vacuum, while for 
the semi-adiabatic system (without geometrical phase), 
mode b is roughly empty. At this instant, the non-empty 
mode is in a Schrodinger cat state. For r = 4A7i n the 
full system has revived; the Q-functions have evolved into 
approximate replicas of their initial states. This is true 
up to an overall phase for the semi-adiabatic case, which 
is typical for a half-revival [53| . From this figure we find 
the revival time for a wave packet encircling the CI in 
the Exe model to be 



T rev = 4ATi„ w 87rA 2 . (32) 



For the semi-adiabatic model, exact revivals (in terms 
of restoring also the overall phase) occur at twice this 
time. It should be pointed out that formula (J22) has 
been verified for a large set of different parameters. 

Even though the phase space distribution of a cavity 
mode is in principle measureable [48], the field intensity 
is directly regained from the cavity output field using a 
photon-counter detector. Already Fig. [5] indicates that 
the average number of photons (n^) differ considerably 
between the full model and the semi-adiabatic one. This 
is verified in Fig. [5] showing the time evolution of (n a ) 
and (rib) for both models. Judging from the field intensi- 
ties in this figure, the revival time of the semi-adiabatic 
model seems to be half the one of the full model, but 
here the a mode is indeed not in a coherent state but in 
a Schrodinger cat. 
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FIG. 5: Snapshots of the Husime Q-functions at times r = 
AT ln (a) and (b), r = 2XT in (c) and (d), r = 3AT ln (e) 
and (f) and r = 4ATi n (g) and (h). The left plots are the 
results with geometric phase, using Hamiltonian (|16[) , and the 
non-geometric phase results, obtained by the Hamiltonian in 
Eq. (|31|l . are displayed in the right figures. The Q-function 
of the a-mode is labeled by a in the plots, while b labels the 
second mode Q-function. The dimensionless parameters are 
A = 6 and Q = 0.5. 



V. CONCLUSIONS 

In this paper we have shown how a system of a two- 
level "particle" interacting with the fields of a bimodal 
cavity may fall in the category of JT models. By repre- 
senting the model Hamiltonian in terms of field quadra- 
ture operators, rather than boson ladder operators, we 
identified its APSs and a CI. In this nomeclature, and in 
particular for the multi-particle analogue (Dicke model), 
the JT effect of cavity QED was identified with the 
normal-superradiant phase transition. The system stud- 
ied here was described by the well known E x e Hamilto- 
nian. Knowledge from earlier research on this model, al- 
most exclusively in molecular or chemical and condensed 
matter physics, has been applied on this cavity QED 
counterpart. Our main interest concerned the geomet- 
rical Berry phase reign from encircling the CI. The effect 
of the geometrical phase was studied by comparing phys- 
ical quantities, such as the field phase space distributions 
and the field intensities, obtained from the evolution of 
either the E x e Hamiltonian or the semi-adiabatic one 



in which no geometrical phase occurs. Clear distinctions 
between the two models were found when the system is 
let to evolve over longer time periods. Energy is swapped 
between the two field modes, and this exchange is highly 

affected by the geometrical phase. From our numerical 

(a) 80r . 1 1 




FIG. 6: The average photon numbers (fu) for both modes a 
(black) and b (gray) as a function of scaled time r/Ti n . The 
upper plot presents the results from using the full Hamiltonian 
(I16p . which includes a geometrical phase, while the lower plot 
shows the results from using Hamiltonian (|31[1 . The effect of 
the geometrical phase is remarkably reflected in the two field 
intensities. The dimensionless parameters are the same as in 
Fig. [5] 



results we could as well present analytical expressions for 
the collapse-revival times for a wave packet encircling the 
CI in the E x e model. 

In addition, by introducing the notion of a wave packet 
evolving on two coupled potential surfaces, a deeper un- 
derstanding of cavity QED models is obtained. This 
work, analyzing the geometrical phase, serves as an al- 
ternative viewpoint of the phenomenon in comparison to 
previous studies such as Refs. [2l|, |22|. It is indeed be- 
lieved that the wave packet method used here will give 
even more thorough insight into cavity QED problems, 
or even trapped ion systems where related CI models are 
expect to occur 55] . We plan to study, using the current 
approach, the dynamics of the Dicke normal-superradiant 
phase transition. Another project underway is to investi- 
gate the " molecular Aharanov-Bohm effect" and "molec- 
ular gauge theory" in terms of cavity QED models. 
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